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PREFACE 


Authority  for  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  to  conduct  this  study.  Work  Unit  No.  31093,  "Harbor  Design  Criteria 
for  Tsunami  Protection,"  under  the  Corps  of  Engineers  Research  and 
Development  Program  was  contained  in  a letter  from  the  Office,  Chief  of 
Engineers  (OCE),  U.  S.  Army,  dated  IT  September  1973.  Funds  were  granted 
by  the  Coastal  Engineering  Research  Area  under  the  field  managership  of 
the  Coastal  Engineering  Research  Center  and  OCE  Technical  Monitor, 

Mr.  Neill  Parker,  HQDA  (DAEN-CWE-H) . 

The  study  was  conducted  by  personnel  of  the  Hydraulics  Laboratory, 
WES,  under  the  direction  of  Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics 
Laboratory,  Dr.  R.  W.  Whalin,  Chief  of  the  Wave  Dynamics  Division,  and 
Mr.  D.  D.  Davidson,  Chief  of  the  Wave  Research  Branch.  The  investigation 
was  conducted  by  Messrs.  A.  W.  Garcia  and  H.  L.  Butler.  This  report  was 
prepared  by  Mr.  Garcia. 

Most  of  the  numerical  computations  were  performed  on  a CDC-7600 
computer  at  the  Los  Alamos  Scientific  Laboratory,  Los  Alamos,  New  Mexico. 

Directors  of  WES  during  the  investigation  and  the  preparation  and 
publication  of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L. 
Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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EFFECT  OF  SOURCE  ORIENTATION  AND  LOCATION  IN  TK 
PERU-CHILE  TRENCH  ON  TSUNAMI  AMPLITUDE  ALONG  Ttt 
PACIFIC  COAST  OF  THE  CONTINENTAL  UNITED  STATES 


PART  I:  INTRODUCTION 


1.  The  U.  S.  Army  Corps  of  Engineers  has  as  an  objective  of  its 
research  and  development  program  the  determination  of  better  harbor 
design  criteria  for  tsunami  protection.  A previous  report1  discussed 
the  subject  of  tsunami  vulnerability  of  the  Pacific  coast  of  the  con- 
tinental United  States  to  tsunamis  originating  in  the  Aleutian  Trench. 
That  report  determined  the  variation  in  tsunami  amplitude  as  a function 
of  coastal  distance  due  to  a standard  uplift  source  at  different  loca- 
tions in  the  Aleutian  Trench.  This  report  is  ^ continuation  of  that 
study  and  addresses  the  subject  of  tsunami  vulnerability  along  the  same 
stretch  of  coast  to  tsunamis  originating  in  the  Peru-Chile  Trench.  In 
addition,  modifications  to  the  numerical  code  used  in  the  previous  report 
allowed  the  simulation  of  the  Chilean  tsunami  of  22  May  i960. 


Historical  Narrat ive 


2.  Five  well-documented  major  tsunamis  have  struck  the  Pacific 

coast  of  the  United  States  in  recent  times:  the  Great  Aleutian  tsunami 

of  19*»6,  the  Kamchatka  tsunami  of  1952,  the  Alaskan  tsunami  of  1957,  the 
Chilean  tsunami  of  i960,  and  the  Alaskan  tsunami  of  196U.  The  196*. 
Alaskan  tsunami  claimed  107  lives  in  Alaska,  b in  Oregon,  and  11  in 
Crescent  City,  California.  In  addition  to  330  lives  in  Chile,  the  i960 
Chilean  tsunami  claimed  6l  in  Hawaii  and  199  in  Japan. 

3.  The  association  of  large  submarine  earthquakes  with  island 

p.L 

arcs  and  deep  submarine  trenches  is  well  documented.  The  Japan- 
Kurile-Kamchatka,  the  Aleutian,  and  the  Peru-Chile  Trench  regions  have 
long  recorded  histories  of  tsunamigenic  activity.  Of  these,  however. 


only  tsunamis  originating  in  the  Aleutian  and  Peru-Chile  Trenches  have 
a history  of  causing  significant  runup  along  the  Pacific  coast  of  the 
United  States.  Because  tsunamis  originating  in  the  Aleutian  Trench  is 
the  subject  of  a previous  report,1  it  will  not  be  recounted  here. 

1*.  On  9 July  1586  a tsunami  occurred  in  the  vicinity  of  Lima,  Peru. 
Waves  as  high  as  81*  ft*  and  inundation  as  far  inland  as  6 miles  were  re- 
ported. On  28  October  171*6  the  city  of  Callao,  Peru,  was  destroyed  by 
two  waves;  reportedly,  only  200  of  5000  inhabitants  survived.  On  20 
February  1835  a major  quake  was  felt  throughout  Chile.  A 50-ft  wave  was 
reported  at  Talcahuano  and  ships  were  sunk  in  the  port  of  Constitucion. 
Waves  were  reported  as  feu*  away  as  Hawaii.  On  13  August  1868  a large 
shock  was  reported  centered  near  Arica,  Peru.  The  city  of  Iquique  was 
reported  inundated  and  waves  rising  30  ft  above  high  tide  were  reported 
at  Talcahuano.  On  11  November  1922  a major  quake  was  reported  near 
Coquimbo  where  three  waves  reached  l-l/l*  miles  inland.  Destructive  waves 
were  reported  from  Japan,  Formosa,  the  Philippines,  and  Hawaii.  On  19 
January  1958  a large  earthquake  of  magnitude  7.5  was  reported  centered 
near  l-l/2°N,  79-l/2°W.  Tsunamis  causing  damage  were  reported  at  Las 
Esmeraldas  and  Guayaquil,  Ecuador,  with  20  killed  and  many  injured."’  The 
tsunami  of  22  May  i960  is  described  as  one  of  the  largest  of  the  past 
century.  Table  1 (adapted  from  Symons  and  Zetler^  shows  a comparison  of 
maximum  wave  heights  for  the  tsunamis  of  191*6,  1952,  1957,  I960,  and 
1961*.  The  comparative  intensity  of  the  i960  tsunami  is  evident,  espe- 
cially when  distance  effects  are  considered;  indeed,  the  generating 
earthquake  was  of  magnitude  8.6,  one  of  the  largest  ever  recorded. 

5.  There  is  some  evidence,  however,  to  indicate  that  Richter 

scale  magnitudes  tend  to  underestimate  source  parameters  (fault  length, 

width,  and  uplift)  and  thus  the  tsunamigenic  potential  of  very  large 

earthquakes.  Cumulative  frequency  versus  surface  wave  magnitude  (M  )** 

s 


A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  and  metric  (SI)  units  to  U.  S.  customary  is 
presented  on  page  3. 

For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 


......  »» »n»i  .» 


plots  show  a linear  trend  for  earthquakes  with  M values  less  than  about 

s 

7-0;  but  for  higher  values  of  M the  locus  becomes  increasingly  steep, 

tending  toward  vertical  at  M approximately  equal  to  8.6.  Therefore,  it 

s 

has  been  commonly  assumed  that  earthquakes  with  M values  greater  them 

s 

about  8.6  do  not  occur,  which  may  indeed  be  true.  This,  however,  does 
not  necessarily  imply  there  is  an  upper  limit  on  earthquake  source  param- 
eters. Conventional  analyses  of  earthquake  energy  spectra  to  determine 
Ms  values  are  frequency-dependent.  Chinnery  and  North7  hypothesize  that 
because  very  large  earthquakes  produce  spectra  which  have  significant 
proportions  of  energy  at  frequencies  below  the  band  where  Mg  is  normally 
measured  and  have  dislocations  which  propagate  for  times  significantly 

greater  than  20  sec,  an  upper  limit  to  M will  be  observed  no  matter  how 

s 

large  the  true  earthquake  size. 
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PART  II:  PERU-CHILE  TRENCH  AS  A TSUNAMIGENIC  REGION 

6.  In  terms  of  plate  tectonics,  the  shallow  seismic  zone  of 

western  South  America  is  a continuous  boundary  along  which  one  plate 

(the  Nazca  plate)  underthrusts  the  adjacent  plate  (the  Americas  plate). 

2 

Kelleher  cites  numerous  studies  that  indicate  this  plate  motion  is 

accomplished  by  large  magnitude  earthquakes;  consequently,  any  segment 

of  the  boundary  that  has  not  ruptured  for  many  decades  should  be  con- 

3 

sidered  a zone  of  high  risk.  Moreover,  Kelleher  has  observed  that 
about  once  each  century  the  entire  fault  segment  near  the  Central  Valley 
province  of  central  and  southern  Chile  (approximately  32°  to  46°S)  has 
fractured  in  a general  north-south  direction  by  a progression  of  large 
earthquakes  and  hypothesizes  that  a new  series  of  quakes  might  start 
about  the  end  of  this  century  near  Valparaiso  (33°S)  and  progress 
southward. 

Q 

7.  Hayes  states  that  crustal  thinning  beneath  the  trench  is 
evidence  that  the  major  tectonic  forces  involved  are  not  laterally  com- 
pressive, and  the  ai  senee  of  these  compressive  forces  can  be  produced 
by  a combination  of  faulting  and  flexure  in  the  transition  region  of 

an  ocean-continent  margin.  Hayes  also  observes  that  almost  all  quake 
epicenters  in  the  area  lie  shoreward  of  the  axis  of  the  trench. 

8.  Records  kept  by  the  Japanese  suggest  that  the  subsea  tectonic 

displacement  of  the  i960  Chilean  earthquake  is  probably  greater  than 

that  associated  with  any  of  the  earlier  earthquakes  during  a period  of 
9 

nearly  500  years.  Indeed,  vertical  movements  of  a larger  regional 
scale  appear  to  be  characteristic  of  that  portion  of  the  trench  extend- 
ing from  Colombia  to  southern  Chile.  The  great  earthquake  of  1835  (para- 
graph 1+)  is  reported  to  have  caused  uplifts  of  as  much  as  8.6  ft  at 
coastal  locations  extending  from  Isla  Mocha  to  Concepcion. 
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PART  III:  TSUNAMI  SOURCE  CHARACTERISTICS 
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9-  Figure  1 shows  the  contours  of  elevation  of  the  initial  sea 
surface  deformation  used  as  input  to  the  simulation  of  the  Chilean 
tsunami  of  22  May  i960.  The  shape  of  the  deformation  is  inferred  from 

Q TO 

the  works  of  Plafker  and  Savage  and  Hwang  et  al.  Plafker  and  Savage 
show  the  maximum  upward  vertical  displacement  of  their  preferred  inter- 
pretation of  the  fault  model  to  be  about  10  m (32.8  ft)  although  the 
maximum  observed  upward  vertical  displacement  was  about  6m  ( 19-7  ft). 

To  avoid  distorting  the  large-scale  features,  the  available  data  were 
smoothed,  i.e.,  the  maximum  vertical  uplift  of  2k  ft  was  chosen  to 
compromise  between  the  maximum  predicted  uplift  of  10  m (32.8  ft)  and 
maximum  observed  uplift  of  6 m (19.7  ft).  Because  the  data  indicate 
there  was  little  downthrow  (less  than  2 m (6.6  ft))  associated  with  the 
tectonic  deformation,  this  feature  was  neglected  in  the  model.  The 
length  of  the  deformation  was  determined  from  the  spatial  distribution 

of  the  epicenters  of  the  initial  quake  and  aftershocks  as  shown  in 
9 

Plafker  and  Savage,  and  both  its  meridional  and  latitudinal  extent  are 

• . 3 

consistent  with  the  area  of  the  rupture  zone  as  shown  in  Figure  2. 

10.  There  is  little  question  that  the  deformation  used  for  the 
simulation  of  the  i960  Chilean  tsunami  is  not  known  with  the  confidence 
of  that  used  to  simulate  the  Alaskan  tsunami  of  28  March  1961+ . Nonethe- 
less, enough  data  are  available  to  credibly  reconstruct  the  bottom  de- 
formation which  generated  the  tsunami. 

11.  Figure  3 shows  the  idealized  axis  of  the  Peru-Chile  Trench 
as  divided  into  12  segments.  A hypothetical  sea  surface  deformation 
which  generates  a tsunami  with  intensity  approximately  equal  to  four  is 
centered  in  and  aligned  with  segments  1,  8,  and  12.  In  order  to  facili- 
tate comparison  with  the  results  of  a previous  study, ^ the  present  hypo- 
thetical uplift  is  identical  with  that  of  the  referenced  study. 
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PART  IV:  C QMFUTAT I ON AL  MODEL 

Equations  of  Motion 

12.  The  numerical  code  used  to  simulate  the  i960  Chilean  tsunami 
and  the  idealized  hypothetical  disturbances  is  based  upon  the  linearized 
long-wave  equations  obtained  by  vertically  integrating  the  Navier-Stokes 
equations  and  equation  of  continuity.  The  equations  as  used  in  the  code 
are  formulated  in  a spherical  system  as  follows: 


au  = in 

3t  R 30 
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(1) 


3V 


in 


3t  Rg  sin  <f>  341 


(2) 


l I 


in  _ _ 

3t  R 


-1 


sin  <j> 


[h {n 


+ d)U  sin  0 + — (n  + d)V 
d<p 


'] 


where 


(3) 


U - depth-averaged  wave  velocity  component  in  the  6 direction 
t = time 

g = acceleration  due  to  gravity 

R = radius  of  the  earth 
e 

n = wave  elevation  from  reference  water  level 
0 = latitude  measured  positive  from  north  pole 
V = depth-averaged  wave  velocity  component  in  the  <j>  direction 
4>  = longitude  measured  positive  eastward  from  Greenwich 
d = local  water  depth 

The  validity  of  using  linearized  equations  and  the  conditions  under 
which  they  hold  are  discussed  in  Reference  1. 


Finite  Difference  Computational  Scheme 


13.  In  order  to  facilitate  simulation  of  the  tsunami  of 
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22  May  i960,  modifications  were  made  to  the  code  used  for  previous 
reports'*" to  permit  larger  grids  and  faster  execution.  In  the 
above-mentioned  reports  the  finite  difference  form*;  of  Equations  1,  2, 
and  3 were  solved  by  an  alternating  direction  multistep  technique  where- 
by during  the  first  half-time  step,  n and  U were  computed  implicitly 
along  lines  of  constant  longitude  and  V explicitly  along  lines  of 
constant  latitude;  in  the  second  half-time  step,  n and  V were  com- 
puted implicitly  along  line  of  constant  latitude  and  U explicitly  along 
lines  of  constant  longitude.  In  the  present  report.  Equations  l*-6  are 
the  finite  difference  forms  of  Equations  1-3. 
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- (Ud  sin  0) 
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(where  it  is  assumed  n<<;d) 

Figure  h shows  the  grid  scheme  used.  The  depth,  d , is  considered 
the  averaged  depth  over  the  grid  square  j to  j + 1 , i to  i+1. 
Equations  H and  5 are  solved  explicitly  at  the  same  half-time  step 
along  lines  of  constant  longitude  and  latitude,  respectively.  Equa- 
tion 6 is  then  solved  for  n one  half-time  step  later.  The  boundary 
conditions,  reflective  at  land  boundaries  and  transmissive  at  the  open 
ocean  grid  boundaries,  are  treated  in  the  same  manner  as  in  Reference  1. 

lH.  The  explicit  code  was  verified  by  comparing  the  results  of 
hindcasting  the  Alaskan  tsunami  of  28  March  196**  with  results  of  hind- 
casting the  same  tsunami  using  the  earlier  implicit-explicit  code. 
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There  was  almost  exact  agreement  of  the  results. 

15-  To  further  reduce  the  computational  time  required  and  to 
allow  a greater  area  to  be  covered,  a 1/3°  by  1/3°  grid  was  used  instead 
of  the  1/3°  by  1/5°  grid  used  in  Reference  1.  The  validity  of  using  the 
coarser  grid  to  propagate  tsunamis  over  the  open  ocean  has  been  demon- 
strated.^ Use  of  the  latest  bathymetric  information  for  the  Pacific 
coast  of  the  United  States  allows  significant  improvement  in  the  detail- 
ing of  the  coastline  and  continental  shelf  area. 

16.  The  depths  for  the  open  ocean  portion  of  the  numerical  grid 

were  interpolated  from  data  of  ocean  depths  averaged  over  areas  of  1° 

lU 

squares  of  latitude  and  longitude  available  from  the  National  Oceano- 
graphic Data  Center.  Normally,  the  numerical  code  interpolates  depth 
data  from  one  point  per  square  degree  to  nine  points  per  square  degree. 
However,  for  the  continental  shelf  area,  the  interpolation  process  was 
not  used;  instead,  data  taken  directly  from  bathymetric  charts  in  a form 
compatible  with  the  numerical  code  were  inserted. 

17.  By  the  modifications  made  to  the  finite  difference  form  of 
the  equations  used  in  the  numerical  code  and  the  use  of  coarser  depth 
grid,  the  computational  time  for  similar  problems  is  reduced  by  a factor 
of  approximately  three  as  compared  with  the  previous  code. 
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PART  V:  PROCEDURE 


18.  To  locate  the  hypothetical  uplifts,  the  Peru-Chile  Trench 
was  divided  into  12  equal  segments  as  shown  in  Figure  3.  An  uplift  of 
the  same  configuration  as  used  in  Reference  1 was  centered  in  seg- 
ments 1,  8,  and  12.  Segment  1,  although  closest  to  the  Pacific  coast 
of  the  continental  United  States,  is  least  favorably  oriented  for  pro- 
ducing large  waves  which  propagate  in  that  direction.  Segment  8,  al- 
though more  favorably  oriented  for  production  of  large  waves,  would  be 
the  segment  most  influenced  by  the  "shadowing"  effect  of  the  reach  of 
coast  immediately  northward  of  the  perimeter  of  the  uplift  area.  In 
principle,  an  uplift  centered  in  segment  7 would  experience  more  shadow- 
ing effect;  but  the  uplift  placed  there  would  have  had  a significant  pro- 
portion on  land  and  therefore  produce  a smaller  tsunami.  The  orienta- 
tion of  segment  12  is  very  similar  to  that  of  segment  8,  i.e.,  more 
favorably  oriented  than  segment  1 for  production  of  large  waves  in  the 
direction  of  the  United  States  but  without  the  shadowing  effect  ex- 
perienced by  segment  8.  One  would  therefore  expect  the  overall  ampli- 
tudes of  the  tsunamis  arriving  at  the  Pacific  coast  of  the  United  States 
to  be  least  for  those  produced  by  the  uplift  centered  in  segment  8 (ex- 
cluding segment  7),  with  tsunamis  originating  in  the  remaining  eight 
segments  producing  values  between  those  of  segments  1 and  12. 

19-  To  compare  the  amplitude  of  the  leading  wave  of  the  tsunami 
generated  by  the  three  hypothetical  ground  motions  as  a function  of  dis- 
tance along  the  Pacific  coast  of  the  United  States,  the  amplitudes  were 
normalized  to  a depth  of  600  ft.  This  is  done  by  analytically  solving 
a one-dimensional  standing  wave  equation. ^ Although  other  methods  for 
calculating  changes  in  wave  amplitude  due  to  changes  in  water  depth 
have  subsequently  been  developed  (e.g.,  numerical  algorithms  based  on 
the  Korteweg-deVries  equations),  it  is  unnecessary  to  employ  them  here. 

In  general,  for  a wave,  the  quantity  (d/L)  is  a measure  of  the  dis- 
persion while  (n/d)  is  indicative  of  the  degree  of  nonlinearity.  For 

p 

a tsunami  at  a depth  of  600  ft,  (d/L)  <<  1 and  (n/d)  <<  1 , indicat- 

ing both  dispersive  and  nonlinear  effects  to  be  small.  The  ratio  of  the 


quantities  |_n/L  (WdrJ  is  termed  the  Ursell  number,  U , and  for 
tsunamis  it  can  be  shown  that  U >>  1 , thereby  making  the  use  of 
Korteweg-deVries  or  Boussinesq  equation  unnecessary.  The  standing  wave 
methodology  is  discussed  in  detail  in  Reference  1 but  will  be  recounted 
briefly  here. 

20.  Combining  the  one-dimensional  long-wave  equation  of  motion 
and  the  continuity  equation  (in  rectangular  coordinates)  yields  the  wave 
equation 


Under  the  assumption  of  simple  harmonic  motion 


and  Equation  7 becomes 


Substituting  d(x)  = d x/a  and  k = u a/gd  where  d is  the  water 

a a a 

depth  at  the  gage  location  some  distance  x = a from  the  shoreline  in 
Equation  9 gives. 


The  solution  to  Equation  10  is 


or  with  the  time  dependence  restored 


Because  a tsunami  has  a very  long  wavelength,  it  can  be  assumed  that 


for  x >_  a . Equating  Equations  12  and  13  yields 


J <2k1/2.1/2> 
o 


Therefore , 

J <2t1/2*1/2> 

n(x)  = C -y  ■ -IV-  cos  (u)t  + e)  (15 

J (2k1/2a1/2) 
o 


From  Equation  15 


J (2k1/,2b1/2) 

n(b) ' c c°9  <“t*£) 
o 


and  substituting  Equation  1 6 


J (2k1/2bl/2) 


n(t)  = — rrrur n(a) 

J (2k1/2a1/2) 
o 


Since  n(a)  is  known  from  the  numerical  simulation,  n at  x = b 

be  determined.  The  values  of  a , b , and  d are  determined 

& 

National  Ocean  Survey  charts. 


21.  Because  the  method  described  above  is  one-dimensional,  it 
accounts  for  the  effects  of  shoaling  and  reflection  but  not  refraction 
or  diffraction.  However,  since  a tsunami  wavelength  is  so  great  in 
waters  deeper  than  600  ft,  few  bathymetric  features  are  large  enough  to 
significantly  refract  or  diffract  the  waves  over  relatively  short  dis- 
tance from  the  nearest  gage  location  to  the  600-ft  contour. 
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PART  VI:  RESULTS 


22.  Table  2 shows  a comparison  of  amplitudes  of  the  leading  wave 
of  the  i960  Chilean  tsunami  as  recorded  by  tide  gages  at  selected  ex- 
posed locations  with  results  of  the  numerical  simulation  at  locations 
nearest  the  tide  gages.  To  avoid  the  false  appearance  of  exact  agree- 
ment at  locations  such  as  Port  Hueneme , the  results  of  the  numerical 
simulation  are  rounded  to  the  nearest  hundredth  foot  rather  than  the 
nearest  tenth.  It  is  quite  proper  to  compare  the  leading  wave  of  the 
actual  tsunami  with  the  leading  wave  of  the  simulation  since,  for  large 
tsunamis  with  minimal  frequency  dispersion,  the  leading  wave  will  prac- 
tically always  be  the  largest  wave  in  water  depths  where  the  comparisons 
are  made  (approximately  600  ft  or  more  at  the  edge  of  the  continental 
shelf).  In  cases  where  subsequent  waves  (at  the  shoreline)  are  larger 
than  the  initial  wave,  it  is  almost  always  due  to  shelf  and  bay  resonance 
which  occurs  shoreward  of  the  point  where  the  comparisons  are  being  made. 

23.  Ideally,  the  amplitudes  resulting  from  the  simulation  should 

be  slightly  less  them  those  recorded  at  corresponding  tide  gages  because 
of  shoaling  of  the  wave  into  shallower  gage  locations.  With  the  excep- 
tion of  the  Presidio  gage,  this  is  indeed  the  case.  Two  factors  which 
could  account  for  the  anomaly  at  Presidio  are:  (a)  frictional  dissi- 

pation of  the  actual  tsunami  over  the  broad,  shallow  shelf  seaward  of 
Presidio  which  is  not  considered  in  the  governing  equations  of  the  code, 
and  (b)  the  averaged  discretized  depths  in  the  finite  difference  program 
cover  too  large  an  area  to  represent  the  detail  of  the  shelf  accurately. 

24.  Figures  5-19  show  contours  of  sea  surface  elevation  of  the 
leading  tsunami  waves  generated  by  simulation  of  the  i960  Chilean  quake 
and  by  hypothetical  uplifts  centered  in  segments  1,  8,  and  12  at  times 
specified  in  the  figures.  The  reason  for  the  extensive  damage  caused 
by  the  i960  Chilean  tsunami  in  the  Hawaiian  Islands  and  Japan  is  evident 
as  the  highest  section  of  the  wave  front  is  directed  toward  these 
locations  (Figures  7 And  8).  The  directivity  of  the  three  hypothetical 
uplifts  is  also  evidenced  in  their  respective  figures,  with  uplifts 

in  segments  1 and  12  radiating  the  largest  waves  primarily  northwestward 


...  .... 


and  segment  o radiating  its  largest  primarily  westward.  The  preferential 
radiational  patterns  of  noncircular  uplifts  is  well  documented. 

25.  Figure  20  shows  the  locations  of  the  64  gages  used  to  deter- 
mine the  leading  wave  amplitude  as  a function  of  distance  along  the 
Pacific  coast  of  the  United  States.  Figures  21-23  are  plots  of  the 
amplitude  of  the  leading  wave  of  tsunamis  generated  by  hypothetical  up- 
lifts located  in  segments  1,  8,  and  12,  respectively.  While  tsunamis 
generated  in  segments  1 and  12  generally  produce  similar  responses  along 
the  Pacific  coast,  the  response  to  the  tsunami  generated  in  segment  8 is 
markedly  less.  This  is  probably  due  to  the  shadowing  by  the  northwest- 
ward to  southeastward  oriented  stretch  of  South  American  coast  between 
lat.  15  and  20°  S.  The  generally  low  response  exhibited  along  the  coast 
from  about  San  Diego  to  Los  Angeles  due  to  tsunamis  generated  in  the 
three  segments  is  probably  the  result  of  sheltering  by  the  Channel 
Islands.  The  peaks  shown  in  the  vicinity  of  Port  Hueneme  (1.4  to  3.3  ft) 
agree  fairly  well  with  heights  recorded  there  during  past  tsunamis,  e.g., 
I960  Chilean,  2 ft;  1957  Aleutian,  <3.5  ft.  The  validity  of  the  large 
peaks  shown  in  Figures  21  and  23  just  south  of  Pt.  Sur  is  unknown. 

No  gage  records  exist  near  that  location  with  which  to  compare  results; 
also,  it  is  a relatively  uninhabited  stretch  of  coastline.  The  re- 
sponse of  the  remainder  of  coastline  is  relatively  uniform  with  a maxi- 
mum amplitude  variation  of  about  1.5  ft.  This  is  not  unexpected  due  to 
the  gradual  eastward  turning  of  the  coastline  as  one  proceeds  northward. 


PART  VII : SUMMARY 

26.  The  Peru-Chile  Trench  has  a long  recorded  history  as  an  area 
of  generation  of  tsunamis  that  cause  significant  runup  along  the  Pacific 
coast  of  the  continental  United  States.  An  explicit  numerical  code  has 
been  developed  and  is  demonstrated  to  adequately  hindcast  the  Chilean 
tsunami  of  22  May  i960.  The  code  is  used  to  simulate  the  generation  and 
propagation  of  tsunamis  with  intensities  approximately  equal  to  four 
which  originate  in  the  Peru-Chile  Trench.  The  amplitudes  of  the  tsunamis 
are  normalized  to  600-ft  depth  and  are  plotted  as  a function  of  distance 
along  the  U.  S.  west  coast  from  the  Mexican  to  the  Canadian  Borders, 
thus  allowing  a relative  comparison  of  the  vulnerability  of  various 
locations  along  the  U.  S.  west  coast  to  tsunamis  originating  in  the 
Peru-Chile  Trench. 

27.  Because  of  priorities  within  the  Civil  Works  R&D  program, 
funding  for  Work  Unit  No.  31093  was  terminated  with  completion  of  the 
work  reported  herein.  If  and  when  funding  of  this  Work  Unit  is  resumed, 
the  logical  next  steps  would  consist  of  (a)  designating  the  locations  of 
existing  and  potential  harbor  sites  on  the  Pacific  coast,  (b)  propagating 
tsunami  waves  into  each  harbor  site  from  probable  tsunami  sources  to 
determine  the  potential  for  damage,  and  (c)  exploring  the  alternatives 

of  relocating  proposed  harbor  sites,  providing  adequate  structural  pro- 
tection, or  modifying  the  hydrography  of  the  areas  in  such  manner  that 
tsunami  waves  would  be  attenuated  to  the  point  that  major  damage  would 
be  prevented.  These  additional  steps  are  necessary  to  meet  the  objec- 
tives of  this  Work  Unit. 
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Table  1 

Maximum  Recorded  Rise  or  Fall  of  Wave  Heights 


Station 

Wave 

Height 

, ft,  for  Year 

1?46 

1?52 

1?57 

I960 

1964 

Honolulu,  Hawaii 

4.1 

4.4 

3.2+* 

5.5+* 

2.7 

Sitka,  Alaska 

2.6 

1.5 

2.6 

3.0 

14.3 

Neah  Bay,  Wash. 

1.2 

1.5 

1.0 

2.4 

4.7 

Crescent  City,  Calif. 

5.9 

6.8 

4.3 

10.9 

(Gage  failed) 

San  Francisco,  Calif. 

1.7 

3.5 

1.7 

2.9 

7.4 

Port  Hueneme,  Calif. 

5.5 

4.7 

3.5 

8.8 

— 

Los  Angeles  (Berth  60),  Calif. 

2.5 

2.0 

2.1 

5.0 

3.2 

La  Jolla,  Calif. 

1.4 

0.8 

2.0 

3.3 

2.2 

San  Diego,  Calif. 

1.2 

2.3 

1.5 

4.6 

3.7 

# Measured  rise  or  fall.  No  extrapolation  made  for  height  beyond  gage 
limit. 


Table  2 

Comparison  of  Elevations  of  Maximum  Rise  of  Initial  Wave 


of  i960  Peru-Chile  Tsunami  to  Numerical  Simulation 


Gage  Location 


La  Jolla,  Calif. 

San  Clemente  Island,  Calif. 

Port  Hueneme,  Calif. 

Presidio  (San  Francisco),  Calif 
Crescent  City,  Calif. 

Neah  Bay,  Wash. 


i960  Peru-Chile 
Tsunami . ft 
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Figure  3.  Idealized  axis 
of  Peru-Chile  Trench  show- 
ing 12  segments 
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Figure  k.  Graphic  representa- 
tion of  the  finite  difference 
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Figure  10.  Contours  of  water-surface  elevation  for  the  leading  tsunami 
waves  9 hr  after  ground  uplift  centered  in  segment  1 


Figure  lb.  Contours  of  water-surface  elevation  for  the  leading  tsunami 
waves  9 hr  after  ground  uplift  centered  in  segment  8 
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Figure  IT.  Contours  of  water-surface  elevation  for  the  leading  tsunami 
waves  4 hr  after  ground  uplift  centered  in  segment  12 
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Figure  19*  Contours  of  water-surface  elevation  for  the  leading  tsunami 
waves  12  hr  after  ground  uplift  centered  in  segment  12 
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amplitude  of  leading  tsunami  wave  generated  in  segment  8 versus  distance 
of  the  recording  gage  from  the  U.  S.  - Mexico  border 
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Wave  amplitude  of  leading  tsunami  wave  generated  in  segment  12  versus  distance 
of  the  recording  gage  from  the  U.  S.  - Mexico  border 
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APPENDIX  A:  NOTATION 

Distance  along  x-axis,  ft 
Amplitude  constant 
Distance  along  x-axis,  ft 
Amplitude  constant 
Depth  of  water 

Water  depth  a distance  "a"  from  shore 
Acceleration  due  to  gravity 

J ( ) Zeroth  order  Bessel  function  of  the  first  kind 

° _i 

Variable,  ft 
Wavelength 

Surface  wave  magnitude 
Refers  to  a time,  nAt 
Radius  of  the  earth 
Time 

Depth-averaged  wave  velocity  component  in  the  0 direction 
Ur sell  number 

Depth-averaged  wave  velocity  component  in  the  <j>  direction 
Distance,  ft 

Length  of  a half-time  step 
Phase  factor 


n Wave  elevation  from  reference  water  level 
0 Latitude  measured  positive  from  the  north  pole 
$ Longitude  measured  positive  eastward  from  Greenwich 

w Wave  frequency,  se 

9 Partial  differential 
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